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■ ABSTRACT 

Aims. The aim of this paper is to study the Gamma Ray Burst (GRB) environment through the analysis of the optical absorption 
features due to the gas surrounding the GRB. 

Methods. To this purpose we analyze high resolution spectroscopic observations (R=20000-45000, corresponding to 14 km/s at 
4200A and 6.6 km/s at 9000A) of the optical afterglow of GRB050730, obtained with UVES@VLT ~ 4 hours after the GRB trigger 
Results. The spectrum shows that the ISM of the GRB host galaxy at z = 3.967 is complex, with at least five components contributing 
to the main absorption system. We detect strong CII*, Sill*, 01* and Fell* fine structure absorption lines associated to the second 
and third component. 

Conclusions. For the first three components we derive information on the relative distance from the site of the GRB explosion. 
Component 1, which has the longest wavelength, highest positive velocity shift, does not present any fine structure nor low ionization 
lines; it only shows very high ionization features, such as C IV and O VI, suggesting that this component is very close to the GRB site. 
From the analysis of low and high ionization lines and fine structure lines, we find evidences that the distance of component 2 from 
the site of the GRB explosion is 10-100 times smaller than that of component 3. We evaluated the mean metallicity of the z=3.967 
system obtaining values x 10^- of the solar metallicity or less. However, this should not be taken as representative of the circumburst 
medium, since the main contribution to the hydrogen column density comes from the outer regions of the galaxy while that of the other 
elements presumably comes from the ISM closer to the GRB site. Furthermore, difficulties in evaluating dust depletion correction can 
modify significantly these values. 

The mean [C/Fe] ratio agrees well with that expected by single star-formation event models. Interestingly the [C/Fe] of component 2 
is smaller than that of component 3, in agreement with GRB dust destruction scenarios, if component 2 is closer than component 3 to 
the GRB site. 
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1. Introduction 

For a few hours after their onset, Gamma Ray Bursts (GRBs) 
are the brightest beacons in the far Universe, offering a superb 
opportunity to investigate both ORB physics and high redshift 
galaxies. Tens of minutes after a GRB, its optical afterglow can 
be as bright as magnitude 13-15. High resolution (a few tens of 
km/s in the optical band), high quality (signal to noise > 10 per 
resolution element) spectra can therefore be gathered, provided 
that the afterglow is observed on such a short time scale by 8m 
class telescopes. 

Fiore et al. (2005) studied the UVES high resolution spec- 
troscopy of GRB020813 and GRB021004 afterglows; they 
showed that the Interstellar Medium (ISM) of GRB host galax- 
ies is complex when resolved down to a width of a few tens of 
km/s, with many components contributing to each main absorp- 
tion system, and spanning a total velocity range of up to thou- 
sands of km/s. Velocity ranges of hundreds of km/s were found 
on the high resolution spectra obtained with MIKE@Magellann 
and HIRES@Keck on GRB051111 and GRB050730 (Chen et 
al. 2005, Prochaska, Chen & Bloom 2006, PCB06 hereafter, 
Penprase et al. 2006). 

The absorption systems can be divided into three broad 
categories. First, those associated with the GRB surrounding 
medium, second, those produced by the ISM of the host galaxy 
along the line of sight but distant enough not to be strongly mod- 
ified by the GRB, and third the intergalactic matter along the fine 
of sight. We discuss these three systems in turn. 

The physical, dynamical and chemical status of the circum- 
burst medium in the star-forming region hosting the GRB pro- 
genitor can be modified by the explosive event, through shock 
waves and ionizing photons. Blueshifted absorbers may be ion- 
ized and even radiatively and/or coUisionally accelerated by the 
GRB and its afterglow. Circumstellar absorption lines in gamma- 
ray burst afterglow spectra can be used to determine the main 
properties of the star progenitors (van Marie, Langer,& Garcia- 
Segura 2005). The variability of the intensity of absorption Unes 
can be used to determine the location and density of the absorber, 
following Perna & Loeb (1998), Boettcher, Fryer & Dermer 
(2002), Draine & Hao (2002), Mirabal et al. (2002), Perna & 
Lazzati (2002), Perna, Lazzati & Fiore (2003). The temperature 
and density of the absorbing gas can be also accurately deter- 
mined by studying particular line ratios, namely those of a fine 
structure transition to its ground state. 

High resolution spectra do not reveal only the circumburst 
matter. They can also be used to probe the ISM of the host 
galaxy along the line of sight. The study of the ISM of z ; 1 
galaxies has so far relied upon Lyman-break galaxies (LBGs) 
at z= 3-4 (see e.g. Steidel et al. 1999) and galaxies which 
happen to be along Unes of sight to bright background quasars 
(QSOs). However, LBGs are characterized by pronounced star- 
formation and their inferred chemical abundances may be related 
to these regions rather than being representative of typical high- 
z galaxies. Faint metal line systems along the line of sight to 
quasars probe mainly galaxy haloes, rather than their bulges or 
discs. Taking advantage of ultra-deep Gemini multi-object spec- 
trograph observations, Savaglio et al. (2004, 2005) studied the 
ISM of a sample of faint K band selected galaxies at 1.4<z<2.0, 
finding Mgll and Fell abundances much higher than in QSO sys- 
tems and sinoilar to those in GRB host galaxies. Such studies wiU 
hardly be extended to higher redshift with the present generation 

* Based on observations collected at the European Southern 
Observatory, ESO, the VLT/Kueyen telescope, Paranal, Chile, in the 
framework of programs 075.A-0603 



of 8m class telescopes, because of the faintness of high-z galax- 
ies. On the other hand, GRB afterglows provide an independent 
tool to study the ISM of high z galaxies. 

Finally, GRB afterglow high resolution spectra can be used 
to probe the Ly-a forest and the high-z intergalactic medium. 
Using GRBs as remote beacons opens up the opportunity to 
highhght any deviation from what is already known from quasar 
forests. For example the so-called "proximity efl'ect" (the feed- 
back of QSOs on their local IGM) should be much reduced for 
GRBs. 

The need to disentangle the contribution to the absorption 
spectra coming from these three different categories renders high 
resolution spectroscopy the ideal technique to study GRB optical 
afterglows. 

Therefore we set up a program to observe bright optical af- 
terglows of promptly localized GRBs with UVES@VLT, taking 
advantage of the capabilities of the GRB-dedicated Swift satel- 
lite. Swift was launched in November 2004, and provides GRB 
positions (2 — 3' precision) in < 10 s. X-ray afterglow posi- 
tions (5" precision) in < 100 s, and an optical finding chart in 
< 300 s (few arcsec to sub-arcsec positional accuracy). These 
characteristics, together with the possibility of quick follow-up 
with UVES@VLT in RRM (Rapid Response Mode), ofi'ers the 
opportunity to obtain high resolution spectra of GRBs with de- 
lays from a few hours down to a few tens of minutes from the 
event. 

To date, five GRBs have been observed with UVES@VLT 
after a Swift trigger: GRB050730, GRB050820A, 
GRB050922C, GRB060418 and GRB060607A. 

In this paper, we present high resolution spectroscopy of 
GRB050730. This GRB was detected by Swift on 2005 July 30 
19:58:23 UT (Holland et al. 2005, Sota et al. 2005) at RA(J2000) 
= 14''08'"17M4andDEC(J2000) = -03''46'l7".8. 133 seconds 
later. Swift' % XRT started to observe the X-ray afterglow in WT 
mode (Barthelmy et al. 2005). Chen et al. (2005) found the red- 
shift of the GRB at z = 3.967 with the MIKE echelle spectro- 
graph on Magellan II, starting the observation ~ 4 hours after 
the burst. The redshift was confirmed by Rol et al. (2005) using 
the ISIS spectrograph on the WilUam Herschel Telescope at the 
Observatorio del Roque de los Muchachos on La Palma. 

Starling et al. (2005) show an early WHT ISIS optical spec- 
troscopy of GRB050730, finding that the host is a low metalUcity 
galaxy, with low dust content. PCB06 present an analysis of the 
MIKE high resolution spectrum of GRB050730, treating the var- 
ious absorption features as single Voigt profiles. Our main goal 
is instead that of trying to disentangle the relative contribution of 
the various components (in velocity space) contributing to each 
absorption system, making a step forward with respect to Chen, 
et al. (2005) and PCB06. 

We focus here on the study of the local medium surround- 
ing the GRB, and in particular on the analysis of the fine struc- 
ture absorption lines produced by the gas in the host galaxy. The 
comparison of such absorption features and their corresponding 
ground states yields information about the gas density, the tem- 
perature and the radiation field. A detailed analysis of the inter- 
vening absorption systems and of the Ly-a forest will be pre- 
sented in a forthcoming paper. 

The paper is organized as follows. Section 2 presents the ob- 
servations and data reduction; Section 3 gives a brief description 
of the absorbing systems identified from the spectra; Section 4 
is devoted to the analysis of the fine structure Unes and to the es- 
timate of the physical quantities characterizing the host galaxy 
gas; in Section 5 the results are discussed and conclusions are 
drawn. 
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Fig.l. The full high resolution spectrum of GRB050730. 
Dashed lines represent the noise spectrum. The spectrum is cut- 
off at the bluemost and redmost end, where the signal-to-noise 
ratio becomes too small. 



2. Observations and data reduction 

In the framework of ESO programs 075.A-0603 we observed 
the afterglow of GRB050730 with the high resolution UV-visual 
echelle spectrograph (UVES, Dekker et al. 2000), mounted at 
the VLT-UT2 telescope. Table 1 gives the log of the observa- 
tions. Both UVES dichroics were used, employing the red as 
well as the blue arm. This allowed us to achieve a particularly 
wide spectral coverage, extending from ~3200A to ~9500A. In 



order to maximize the signal to noise ratio the CCD was re- 
binned by 2 X 2 pixels. The data reduction has been performed 
using the UVES pipeline (Ballester et al. 2000). The final useful 
spectra extend from about 3750A to about 9780A and were re- 
binned to 0.1 A. The resolution element, set to two pixels, ranges 
then from 14 km/s at 4200A to 6.6 km/s at 9000A. The noise 
spectrum, used to determine the errors on the best fit line param- 
eters, has been calculated from the real, background subtracted 
and rebinned spectrum, using line free regions. This therefore 
takes into account both statistical errors and systematic errors in 
the pipeline processing as well as the background subtraction. 



3. Absorbing systems 

Fig. 1 shows the high resolution spectrum of GRB050730, with 
a cut-ofi" at both the bluemost and redmost end, where the signal- 
to-noise ratio becomes too small. The spectrum exhibits a large 
number of absorption features, belonging to different absorption 
systems. The main system is due to the gas in the host galaxy, 
at redshift zcrb - 3.96764. This value confirms the previous 
redshift deternoinations for this GRBs. We identified four more 
intervening absorbers at z < zgrb not related to the GRB host 
galaxy but rather to galaxies along the line of sight. Each system 
is briefly described in the next subsections. 

3. 1. The host environment atz = 3.967 

This is the system with the greatest number of features on the 
GRB line of sight. Fig. 1 clearly shows the Ly-a absorption at 
~ 6000 A. The redshift inferred from this feature is confirmed 
by a large number of metal-line transitions. The precise value, 
which we calculated from the analysis of the narrowest lines 
associated with the host galaxy, is z = 3.96764. The spectrum 
exhibits a large number of absorption fines, from neutral hydro- 
gen (Ly-a, Ly-yS, Ly-y), via neutral metal-absorption lines (O I) 
and low ionization lines (Cll, Sill, AlII, Nill, Fell, PII, SII) 
to high-ionization absorption features (CIV, SiIV,N V, O VI). 
In addition, strong fine structure fines (CIP, Sill*, OI* Oil**, 
Fell*, Fell**, Fell***) have been identified, meaning either 
that the GRB environment is composed of a high density gas or 
that an intense radiation field is exciting such features. 

The ISM of the host galaxy is complex, with many compo- 
nents contributing to the absorption system. This complex envi- 
ronment, together with a detailed study of the fine structure lines 
aimed at estimating the physical quantities of the circumburst 
gas, will be described in the next section. 

3.2. The intervening absorber atz- 3.564 

This system presents neutral hydrogen features (Ly-a, Ly-jS), 
neutral metal-absorption fines (CI, Si I), low ionization fines 
(Al II, Mg II, Si II, Fe II) and high ionization lines (C IV, Si IV, 
N V). The redshift calculated from these features is z ~ 3.56395. 

3.3. The intervening absorber atz = 2.262 

This system presents only a few neutral metal-absorption lines 
(OI, Sil, All) and some low ionization transitions (AlII, Fell). 
The redshift of this system is z ~ 2.26181. 
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Table 1. Journal of observations 



Date UT 


Dichroic 


B. ArmC.W. 


R. Arm C.W. 


slit" 


seeing" exposure* 


time since GRB'' 


07/31/05 00:07:17 


1 


3460A 


5800A 


1 


;i 50 


4.15 


07/31/05 01:02:09 


2 


4370A 


8600A 


1 


;i 50 


5.07 



"arcsec; ''min; '^hr 



3.4. The intervening absorber at z - 2.253 

This is actually a double absorption system, the redshifts of the 
two components being z ~ 2.2526 and z ~ 2.2536. It shows neu- 
tral metal-absorption lines (Si I, Mg I), low ionization transitions 
(Si II, Fe II, Mg II) and possibly a C IV high ionization line. 

3.5. Tlie intervening absorber at z - 1.772 

This is another double absorption system, the redshifts of the 
two components being z ~ 1.7723 and z ~ 1.7729. The system 
shows neutral metal-absorption lines (Mg I, Fe I) and low ion- 
ization transitions (Mg II, Fe II). 

4. The main system at z=3.967: physical properties 
of the GRB environment 

This section, as well as the rest of the paper, is devoted to the 
study of the main absorption feature of the spectrum, produced 
by the gas in the host environment, located close to the GRB 
explosion site. 

The analysis of the GRB environment is intricate due to the 
complexity of the absorption lines of the spectrum, which in sev- 
eral cases can hardly be fitted with a single Voigt profile. This 
means that many components contribute to the gas in the GRB 
environment. In other words, several shells of gas which may 
be close to or far from each other, appear mixed together in the 
spectrum (in velocity space). Our main goal is to disentangle 
their relative contribution. 

The identification of the components of the absorption lines 
and the fitting procedures will be discussed in the next subsec- 
tion. 

4. 1. Absorption features of tlie GRB environment 

The gas surrounding the GRB constitutes a complex environ- 
ment, and this is reflected in the complexity of the absorption 
features. The presence of several components is indicative of 
clumpy gas in the GRB environment, composed of different ab- 
sorbing regions each with diff'erent physical properties. For ex- 
ample, a first component is present for the high ionization lines 
O VI and (less intensely) C IV, but not for Si IV and lower ion- 
ization lines (Fig. 2). This means that this particular component 
probably comes from a region very close to the GRB site, where 
it experienced a very intense flux of radiation ionizing the atoms 
several times. 

If the absorption features are treated as a single Voigt profile, 
we obtain very bad fitting parameters and very high values. 
The identification of the different components is somewhat sub- 
jective, the true message being that the geometry and kinematics 
of the ISM clouds probed by the GRB line of sight are complex. 
We decided to use the C IV /11548,1550 A and Si IV /11393,1402 
A doublets, which have the wider velocity range, to identify the 
different components constituting the circumburst matter. Fig. 2 
shows CIV /11550 A at z = 3.9676. For this line we obtain a 
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Fig. 2. The C IV /11550 A and Si IV ^1402 A absorption features 
together with the O VI /i 103 1,1 037 A doublet. Each component 
is identified with progressive numbers from the highest wave- 
length to with the lowest one. 



good fit (x^ - 1 .04) if we assume that five different components 
contribute to the absorption features. 

In the following, the components are numbered so that the 
first component is the one with the longest wavelength and the 
highest, positive, velocity shift with respect to a zero point ar- 
bitrarily placed at z=3. 96764, and so on. We next fixed the red- 
shifts of the C IV components and tried to fit the Si IV doublet 
with the same five Voigt profiles. We did not obtain a good fit for 
the following reasons. First, the Si IV doublet does not exhibit 
the first and the last components that we see in the CIV (i.e., 
there are only three components), and second, a single central 
component is not enough to model the central part of the Si IV 
features alone. More precisely, the central component (number 
3) must be splitted in two parts (3a and 3b) in order to obtain a 
good fit for Si IV; infact, such a feature does not show satura- 
tion in its central part, while CIV does. This is clearly visible 
especially for Si IV /11402 A (Fig. 2). Thus, using four compo- 
nents to represent Si IV, the first and the last being fixed at the 
redhift of the second and the fourth of C IV, we obtain a good fit 
(X^ = 0.64) for this feature. All the further features we analyzed 
have been fitted fixing the redshift of the components to those we 
identified for either C IV or Si IV. Case by case, we decided to fit 
groups of elements with either C IV or Si IV, the choice depend- 
ing on the extension of the lines and the saturation of the central 
part of the features. We specify that in these fitting procedures, 
the redshift of the components has been linked to that of C IV or 
Si IV, which in turn has been fixed to the values we find in the 
C IV and Si IV fit. Similarly, the column densities and Doppler 
parameters of the components of C IV and Si IV has been fixed 
while fitting these pilot lines together with other ones. 

Finally, since the two central components of Si IV (3a and 
3b) can not be resolved in the case of C IV, the corresponding 
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Fig. 3. The C II /11036,1037* A fine structure doublet. The soUd 
line represents the five Voigt components fit. 



column densities are computed separately, but summed together 
in the analysis as if we were dealing with a single component. 

We detected several absorption fine structure features in the 
spectrum at the GRB redshift. The gas can be excited to such 
states via collisional effects (if the density is sufficiently high) or 
by the absorption of UV photons followed by the spontaneous 
radiative decay to an excited level of the ground state. Such lines 
can be very important in determining some physical properties of 
the gas, and thus we will consider them carefully. The next sub- 
subsections describe the fits we performed to obtain the column 
densities of the elements that show fine structure transitions. The 
column densities corresponding to each component are reported 
in Table 2 and 3. 

4.1.1. The ion CM 

For this ion, we identified two fine structure doublets, namely 
CII /11036 A CII* ^1037 A and C II ^1335 A CII* ^1337 A. 
We fitted these doublets together with CIV (Figs. 3 and 4). We 
note that only the second and third components out of the five 
in the fits in Figure 2 are present in the excited fine structure 
levels, with a possible indication of the first component for C II* 
/1 1037 A. The fifth component, which is the one with the lowest 
redshift, is present only for C II /1 1036 A; the blending of C II* 
/1 1037 A by C II /1 1036 A. makes impossible to say whether this 
component is present in the excited level or not. The C II /11335 
A doublet appears to be saturated; this is evident by looking at 
the best fit, which is completely flat in the center of this line 
complex. Therefore, we use only the CII /11036 A doublet to 
evaluate the electron density and temperature. 

4.1.2. The ion Sill 

For this ion, we identified three fine structure doublets: Sill 
/11304 A (Sill* ^1309 A), Sill /11526 A (Sill* ^1533 A) and 
Sill /11260 A (Sill* ^1265 A). We fitted the first two doublets 
together with Si IV (Figs. 5 and 6). The first of these doublets 
has been fitted with the OI fine structure triplet too. 

The third doublet falls in the Ly-a absorption gap. A reli- 
able analysis of this feature is difficult because a correct sub- 
traction of the continuum is not possible. Moreover, the excited 
level is further splitted into two sublevels, which blend with each 




-ISO -BO BO 160 

kJelDcity CKm/s) 



Fig. 4. The CII ^1334,1335* A fine structure doublet. The solid 
line represents the five Voigt components fit. 
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Fig.5. The Sill ^1304,1309* A fine structure doublet and OI 
/il302,1304*,1306** A fine structure triplet. The solid line is 
the four Voigt components fit. 



other and render it impossible to fit the contribution of the vari- 
ous components. 

As for CII, only the second and possibly the third compo- 
nent is present in the excited fine structure levels. Finally, the 
third component of the Si II* /11533 A, if present, would fall in a 
region rich of telluric features, and thus can hardly be useful in 
determining the physical parameters of the gas. 

4.1.3. The atom 01 

This atom is present with the OI ^1302,1304*, 1306** A fine 
structure triplet (Fig. 5). It has been fitted with a four compo- 
nent profile, together with the Si IV components and the Si II 
/II 304, 1309* A fine structure doublet. As for previous ions, the 
only components present in the excited fine structure levels are 
the second and the third (splitted in 3a and 3b). The second com- 
ponent appears to be saturated both in the ground and in the first 
excited level O I /11304* A, and will not be used in the analysis. 
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Fig. 6. The Si II /11526,1533* A fine structure doublet. The solid 
line is the four Voigt components fit. The /11533 line is blended 
by several telluric features. 
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Fig. 7. The Fell 1608 fine structure multiplet, together with the 
SilV 1393 absorption feature. The solid line represents the dou- 
ble Voigt component fit. 



4.1 .4. The ion Fell 

This ion is present in a wide variety of flavors. In this paper, 
we concentrate on the fine structure multiplet around A - 1600 
A. This multiplet shows the following absorption features: Fe II 
/11608 A,FeIU1611 A,FeII*/ll618A,FeII*.11621 A,FeII** 
^1629 A, Fell*** il634 A, Fell*** /11636 A. Fig. 7 presents 
the multiplet, together with the Si IV /II 393 A absorption fea- 
ture. 

We can see that the Fe II excited levels only show the second 
component, while the ground state also includes the third com- 
ponent. We obtain a good fit with a single Voigt profile for the 
excited levels and a double line profile for the ground state. The 
central wavelengths of these Voigt profiles have been fixed to the 
redshifts of the second and third components of C IV. 

In the next subsections we will describe in detail the physics 
which relates the fine structure transitions to the density and tem- 
perature of the gas that produces these features (provided that the 
fine levels are excitated by collisional effects). 



4.2. Fine structure line excitation via coliisions 

For ions with a doublet fine structure in the ground state, the 
ratio between the population density of the higher and the lower 
level is: 



"2 _ Q\2 + Wi2 
«1 221 + W21 ' 



(1) 



where Qn and Q21 are the collision rates from level 1 to level 
2 and vice versa, and wn (w2\) is the photon absorption (decay) 
rate between the two levels (Bahcafl & Wolf 1968). Qn and Q21 
can be written in the form Q - < crv >j Nj, where each 
kind of exciting particle (proton, electron, etc.) is denoted by j, 
cr is the collisional cross section for the process involving the j- 
th species, whose density and average speed are Nj and < v >y, 
respectively. 

We are interested in the situation in which level 2 is the ex- 
cited fine structure of the ground state 1 . The excited fine struc- 
ture levels can be produced by photon pumping or collisional 
excitations. Generally, the following mechanisms can be respon- 
sible for this excitation: (i) direct excitation by IR photons; (ii) 
excitation to upper levels through absorption of UV photons, fol- 
lowed by the spontaneous decay to an excited level of the ground 
state (indirect UV pumping); (iii) collisions with charged parti- 
cles (mainly electrons) or with neutral hydrogen. PCB06 showed 
that the excitation rate by UV pumping is larger than IR pump- 
ing, rendering the latter excitation mechanism unimportant for 
gas observed in GRB afterglow spectra. In the following, we will 
elaborate the scenario in which collisions are the main mecha- 
nism responsible for the production of fine structure transitions 
and discuss subsequently (section 4.4) the situation in which in- 
direct UV pumping is at work. 

If we assume that the collisions between electrons and ions 
represent the main process responsible for the excitation of the 
fine structure levels, equation (1) simplifies since the quantities 
of the parameter Q now refer only to electrons. In particular, for 
the C II and Si II fine structure doublets, equation (1) becomes: 



75 



rie - 



[2 



1] 



(2) 



and 



1.28 X 10^^ 



(3) 



respectively (Srianand & Petitjean 2001). These equations have 
been derived taking into account that the radiative excitation 
from the ground to the excited state is forbidden, and using the 
statistical weights; the radiative decay rates and the energy gaps 
are given in Bahcall & Wolf (1968). Here is the electron den- 
sity, is the column density of the element in parentheses (the 
asterisk refers to excited fine structure levels) and T4 is the ki- 
netic temperature in units of 10'* K. 

For a fine structure system with multiple levels, the situation 
is more complicated, unless we assume that electron-ion colli- 
sions are dominant with respect to the radiative processes. If this 
is the case, the w terms in equation (1) can be neglected, and the 
electron density disappears, yielding: 



^ = ^ exp{-EijlkT], 
nj gj 



(4) 
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Fig. 8. Top panel: the electron density as a function of the ratio 
between the absorbing columns of the ground state with respect 
to the excited one, for three different values of the temperature, 
calculated for the C II fine structure doublets. Vertical lines rep- 
resent the observed (N/N*) ratio (solid Une) and its Icr errors 
(dashed lines). The figure refers to the second component of the 
host environment gas. Bottom panel: same as top panel, but for 
Si II fine structure doublets. 



in the case of thermal equilibrium (PCB06). Equation (4) tells 
us that the levels are populated according to a Boltzmann dis- 
tribution; here gi is the degeneracy of the state i, Eij = £, - Ej 
is the energy difference between the state i and j and T is the 
excitation temperature. 



4.3. Temperature and density of the gas in the case of 
collisional excitation 

In this section we try to obtain an estimate of the temperature 
and density of the gas in the GRB environment. We showed 
in the previous section that the excited states of the fine struc- 
ture transitions present some absorption only in the second and 
third component of the five identified. For this reason, the phys- 
ical quantities can be computed or constrained only in these two 
cases. 

We proceed in the following way. First of all, we assume that 
the electron-ion collisions are the main mechanism responsible 
for the production of the excited fine structure feature (PCB06). 

For fine structure doublets (C II and Si II), we plot in Fig. 8 
the electron density as a function of the ratio between the ab- 
sorbing columns of the ground state with respect to the excited 
one (N/N*), using eqs. (2) and (3), and for different values of the 
temperature. In this plot, vertical lines are loci of constant N/N*. 
We then draw vertical lines corresponding to our measure of the 
ratio N/N* and to the corresponding errors. In this way we can 
constrain the density of the gas for a fixed interval of tempera- 
ture, and even constrain the temperature itself, if both the CII 
and Sill present an excited fine structure absorption. 

For fine structure multiplets, we use eq. 4 under the assum- 
pion of thermal equilibrium. Fitting the ratio of the column den- 
sities to their multiplicity factors with respect to the energy dif- 
ference between each excited state and the ground state, we can 
obtain a measure of the temperature of the gas. Such a tempera- 
ture can constrain the gas density using the plots of eqs. (2) and 
(3) for the fine structure doublets, or can be used as an input for 
the PopRatio software package (Silva & Viegas 2001), to esti- 
mate the electron density independently once the ratio between 
the column densities of the fine structure levels of the multiplets 
is known. 

As previously remarked, only component 2 and 3 of the main 
absorption system feature fine structure transition; in the follow- 
ing sub-subsections we present the analysis of such components. 



4.3.1. Component z=3.967_2 

This component features high column densities for all the 
ground states: C II, O I, Si II and Fe II. The excited fine struc- 
ture levels appear in all their flavors: we have strong C IP, O I* 
/ OP* and Si IP columns; moreover Fell is present with its 
first three excited levels. We plot the electron density as a func- 
tion of (N/N*) for three different temperatures, namely T = 10'^, 
r = 5 X 10^ and r = 10^ K, for C n (Fig. 8, left panel) and Sill 
(Fig. 8, right panel). For CH we used the 1036-1037 A feature, 
because C II /11334 A is saturated. 

We can see from Fig. 8 that both for C II and Si II the ob- 
served column densities can only put a lower limit to the elec- 
tron density, independently of the temperature of the gas. This 
lower limit lies in the range 50 200 cm"^ and 100 h- 300 cm"^ 
for temperatures of 10^ lO"' K, for Cn and Sin respectively. 

For this component, we also made a fit to the temperature 
using eq. (4) and the Fell multiplet data (Fig. 9). The value we 
found is r = 26001^000 k (90% confidence limit); such a value 
agrees with that found by PCB06, for the same GRB. 



The same fit has not been performed for the O I, because its 
ground state appears to be saturated. 
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Table 2. Logarithmic ion column densities in cm~^ for the components of the main absorption system at z=3.967 



Velocity shil'l" (coniponcnl) 


Cll 


CIV 


01 


Sill 


SilV 


Ic 11 


+32.6 (1) 


Blend 


14.4 ± 0.08 


NO 


NO 


NO 


NO 


+2.4 (2) 


14.33 ± 0.15 


14.86 + 0.08 


SAT 


14.68 + 0.23 


15.53 + 0.11 


15.24 + 0.06 


-44.0 (3a) 


14.81 + 0.13 


14.75 + 0.12 


14.94 ± 0.09 


13.62 + 0.06 


13.72 + 0.09 


14.66 + 0.22 


-44.0b (3b) 






13.43 + 0.11 


13.60 + 0.07 


14.11 +0.10 




-90.2 (4) 


12.98 ± 0.3 


13.57 + 0.13 


NO 


12.53 + 0.24 


13.21 + 0.03 


NO 


-154.6 (5) 


12.75 ± 0.35 


13.04 + 0.06 


NO. 


NO 


NO 


NO 



'km/s 



Table 3. Logarithmic ion column densities of excited fine structure states in cm ^ for the components of the main absorption 
system at z=3.967 



Vel." (Coinp) 


C 11* 


Ol* 


O 1** 


Si 11* 


I'C 11* 


l-cll** 


I'C 11*** 


+2.4 (2) 


14.69 + 0.26 


SAT 


14.64 + 0.05 


14.44 + 0.09 


14.29 + 0.03 


13.74 + 0.08 


13.60 + 0.09 


-44.0 (3a) 


14.55 + 0.09 


13.56 + 0.08 


13.46 + 0.09 


13.41 +0.15 


< 13.14 


NO 


NO 


-44.0b (3b) 
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12.55 + 0.60 


13.66 + 0.11 


NO 
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'km/s 
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Fig. 9. The ratios between column density and degeneracy of the 
excited fine structure Fe II levels as a function of their energy 
gap to the ground state. The solid line represents our best fit to 
the data. 

4.3.2. Component z=3.967_3 

This component features high column densities for the ground 
states of the C II, O I and Fe II, while that of Si II is very low. 
As for the excited fine structure levels, we have strong C II* and 
low OI columns; Fell excited levels are absent. Sill* is dilR- 
cult to treat. Apart from the -11265 line that falls in the Ly-a 
absorption gap, the /11533 line is blended with several telluric 
features (Fig. 6). The last flavor of Si II, the /II 309 line, shows a 
hint of a possible absorption feature (Fig. 5); however, since the 
corresponding ground state is very lowly populated, this feature 
is somewhat suspect. As for the second component, we plot the 
electron density as a function of (N/N*) for three different tem- 
peratures, namely T = 10^ J = 5 x 10^ and J = 10^ K, for C II 
and Sill (Fig. 10). For CII we used the A 1036-1037 feature, 
because C II /I1334 A is saturated. 

Fig. 10, left panel indicates that the electron density inferred 
from CII lies in the range 10 ^ 60 cm"^ for a temperature of 
10^ -H 10"* K. If we consider the plot for Sill, the density appears 
to be higher than ~ 200 cm"^, if we assume a temperature of 
~ 10^ K. Thus, there is no overlap between the two solutions. 
As we discussed above, the Sill* /11309 A is a suspicious fea- 



ture, because the corresponding ground state (Sill AI304 A) is 
poorly populated, while for the other elements the ground state 
is very highly populated when an excited level is present. In the 
following, we will consider either that Sill* /11309 A is real or 
it is not. 

For this component, we also estimate a temperature using eq. 
4 and OI triplet data. The value we find is T - 10^ K, but we 
only have two points in the fit so the errors are too large to put 
an acceptable constraint on the temperature. 

4.4. Fine structure line excitation via indirect UV pumping 

In the previous sections we assumed that fine structure excited 
levels are produced by collisions. The other, competitive mech- 
anism is indirect UV pumping. In fact, the fine structure levels 
can be excited by the absorption of a UV photon to an upper level 
followed by the spontaneous decay to an excited lower level. If 
this mechanism is at work, we cannot use the fine structure abun- 
dances to gather information about the temperature and density 
of the gas, but only on the strength of the radiation field. PCB06 
(using the PopRatio code developed by Silva & Viegas 2001, 
2002) analyzed the relation between the far UV radiation field 
intensity and the relative fraction of excited fine structure states 
with respect to their ground levels (for OI, Sill and Fell, see 
their Figs. 7 and 8). In the following subsections, we will use 
their plots to estimate the radiation field intensity from our col- 
umn density data, in the case that the fine structure transitions 
are produced by indirect UV pumping. 

4.4.1. Component z=3.967_2 

As we showed before, this is the component which presents more 
fine structure excited states. We evaluate the ratio between the 
ground and the first excited levels for FeII(/ll608-1611 Aand 
AI6I8*-I62I* A) as well as flie ratio between the Sill /II526 
A and Sill* /11533 A fine structure levels. The OI ground state 
is saturated and can not be used in the analysis. Using Fig. 7 of 
PCB06, we find that the two ratios are compatible with a radi- 
ation field intensity of G/Gq = 5 X 10^, Go = 1.6 x 10"^ erg 
cm"^ s ' being the Habing constant. A consistent result is ob- 
tained estimating the ratios of the absorbing column of the first 
excited level of Fell to the second and third ones, and using Fig. 
8 of PCB06. Again, we find that the two ratios yield a value of 
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section 4.3.2, the Sill* detection is uncertain. Nevertheless, us- 
ing the O I and Si II ratios between the ground and first excited 
level, we find two different values for the radiation field inten- 
sity: G/Go - 10^ and G/Gq = 10^, respectively. These values 
do not overlap at the 90% confidence limit. Since the detection 
of Si II* is uncertain, we cannot exclude in this way that the fine 
structure excited levels in this component are produced by in- 
direct UV pumping. However, the theoretical lines representing 
the FeII*/FeII and OI*/OI abundance ratios as a function of 
the radiation field in Fig. 8 of PCB06, approach each other very 
closely. This means that we should observe a similar column of 
Fell* with respect to Fell, as we observe for OI* with respect 
to O I. Since Nqi* /Nqi = 0.04, and Npeii = 14.66 cm'^ (see Tab. 
2), we expect A^fe//* - 13.26 cm"^. However, this is not the case, 
because the upper hmit that we found for the first excited level of 
the Fell is Npeii* < 13.14 cm"^ at the 90% confidence level. On 
the other hand, if the fine structure states of the third component 
are produced by collisions, the O I*/0 1 abundance ratio would 
yield an electron density of a few tens cm"^, compatible with 
that found for the C II fine structure doublet in section 4.3 .2, and 
the ratio FelP/Fell would be 10 times smaller the OP/OI one 
(See Fig 9 in PCB06). This means that we expect = 11.96 

cm"^, a value consistent with the observed upper limit. 



1000 




'""'sill/'^Sill 



Fig. 10. Top panel: the electron density as a function of the ratio 
between the absorbing columns of the ground state with respect 
to the excited one, for three different values of the temperature, 
calculated for the C II fine structure doublets. Vertical lines rep- 
resent the observed (N/N*) ratio (solid line) and its Icr errors 
(dashed lines). The figure refers to the third component of the 
host environment gas. Bottom panel: same as top panel, but for 
Si II fine structure doublets. 



G/Go = 5 X 10^ The 90% confidence interval for the field mten- 
sity is G/Go = 10^ ^ 10^ 



4.4.2. Component z=3.967_3 

This component exhibits the C II fine structure doublet, the O I 
triplet and possibly the Si II excited level. As we discussed in 



4.5. Neutral elements and gas distance from the GRB 

The presence of neutral elements in the spectra of GRB after- 
glows is important, because it can put strong constraints on the 
distance of the gas from the GRB. This is because the first ion- 
ization potential of the elements usually lies below the 1 Ryd 
threshold. Thus, the intense radiation coming from the GRB 
would easily ionize the neutral elements, without the screening 
of the Hydrogen atoms, the first ionization potential of which is 
at higher energies. PCB06 observed the Mg I A2852 A transition 
in GRB051 1 11, and they put a lower limit of 50 pc to the dis- 
tance of the gas from the GRB. Moreover, they refined this limit 
to 80 pc since no variability of the Mg I line has been observed. 

In GRB050730, we observe many O I transitions. Oxygen, 
however, has a first ionization potential that lies just above the 1 
Ryd threshold, and can exist in its neutral state even in the pres- 
ence of a strong radiation field, because the Hydrogen screens 
it from the UV photons. We can not observe the Mg I /12852 A 
transition because of the redshift of our GRB. In addition, we 
do not observe any other neutral element in the afterglow so we 
can not put any constraint on the gas distance from the GRB 
this way. The upper Umit that we find for the neutral carbon ab- 
sorption line (CI /11656 A) is Na* < 13.2 cm"^ at the 90% 
confidence level. 

Some considerations on the relative distance of the difi'erent 
shells of gas from the GRB can be pursued. Component 1 of the 
z = 3.967 system, the one with the longest wavelength, highest 
positive velocity shift (see sec. 4.1), does not present any fine 
structure or low ionization lines; it only shows very high ioniza- 
tion features, such as C IV and O VI, suggesting that this compo- 
nent is very close to the GRB site. Moreover, if we assume that 
the fine structure lines of the shells are produced by the same 
mechanism (collisions or UV pumping), we can also say some- 
thing about component 2 and 3. In fact, in the case of collisional 
excitations we, can estimate the electron density of the gas and 
this can be related to the distance r using the definition of the 
ionization parameter: 
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Fig. 11. The Ly-/3 and Ly-y absorption features, together with the 
C IV /11548 and C IV A1550 doublet. The red Une represents the 
five Voigt component fit. 

where Lp/,o, is the number of ionizing photons and c is the speed 
of light. A time independent estimate of U with photoionization 
codes like CLOUDY (Ferland 2002) is inappropriate in the case 
of a strongly variable phenomenon such as a GRB. Nevertheless, 
if we assume the electron density found in section 4.3 for the 
third component (10 -H 100 cm"^) and for the second one (10^ H- 
10* cm"^) we can conclude (using the scaling law of eq. 5) that 
the gas belonging to the third component lies 10 -H 100 times 
farther from the GRB than that belonging to the second one. 

In a similar way, we can assume that indirect UV pumping is 
at work. In this case, the distance increases as the square root of 
the UV flux. If we trust the estimate of the radiation field inten- 
sity evaluated only from the O I*/0 1 ratio for the third compo- 
nent, we can conclude that the gas in the second component is a 
few ten times closer to the GRB than that of the third one. 

4.6. Metallicity 

Fig. 1 1 shows the Ly-/3 and Ly-y absorption features. It is im- 
possible to discriminate between the five components identified 
from the analysis of C IV and Si IV. The hydrogen absorption 
features result in very broad lines with two extended wings, 
which reach well beyond the velocity range of the heavier ions. 
Since a component by component analysis is not possible, we 
can only estimate the metallicity of the circumburst material in 
terms of the ratio [X/H] between the total absorption column of 
the X element with respect to that of hydrogen. Since the latter 
mainly includes hydrogen in the outer regions of the galaxy, we 
will underestimate the metallicity. 

Table 4 shows our estimate of the metallicity. The second 
column presents the column densities of the different ions, be- 
ing the sum of the columns of the five components of the ground 
states and of the excited fine structure levels, if these are present. 
The third column of Table 4 is the total column density for each 
element, and is obtained adding the column densities of the vari- 
ous ions through which the element appears in the spectrum. We 
could not get an estimate of the column density of Oxygen, be- 
cause the absorption features of O I and O VI are saturated. Not 
surprisingly, we obtain a metallicity of 10"~ lO"-' relative to the 



solar one, in agreement with the paper by Starling et al. (2005). 
It should also be noted that these estimates are not corrected for 
depletion in dust grains. In the Galactic ISM only ~ 1% of the 
iron is in the gas phase while the rest is locked in dust grains 
(Savage & Sembach 1996). The situation can be very different 
in a GRB host galaxy, and in particular, in the GRB surrounding 
medium (Waxman & Draine 2000, Draine & Hao 2002, Perna & 
Lazzati 2002, Perna, Lazzati & Fiore, 2003), because the intense 
GRB radiation field can efficiently destroy dust grains. Caution 
should therefore be used in using mean metallicity estimates. 

We can also calculate the ratio [C/Fe], representative of the 
enrichment of the a elements relative to iron. We measure a 
mean [C/Fe]=0.08±0.24, consistent with the value predicted by 
the models of Pipino & Matteucci (2004, 2006) for a galaxy 
younger than 1 Gyr subject to a burst of star-formation. In such a 
case Pipino & Matteucci (2006) also predict a low [Fe/H] value, 
close to that in Table 4. Again, we warn that mean abundance 
estimates may be affected by large systematic uncertainties. 

A safer approach is to estimate the [C/Fe] of each component 
of the z=3.967 system. We find that [C/Fe] of component 2 is 
-0.15+0.13 while that of component 3 is +0.53+0.23. Indeed, 
note in Tables 2 and 3 that while the total carbon column density 
of component 2 and 3 are similar, the total colum density of iron 
of component 3 is about four times less than that of component 2. 
A similar conclusion applies to silicon. The total silicon column 
density of component 3 is also about 10 times less than that of 
component 2. Interestingly, Perna, Lazzati & Fiore (2003) found 
that silicates tend to be destroyed more efficiently than graphite 
if a dusty medium is exposed to the intense GRB radiation field. 
This would leave more iron free in the gas phase in clouds closer 
to the GRB site than in farther clouds, in agreement with our 
suggestion that component 2 is 10-100 times closer to the GRB 
site than component 3. 

5. Conclusions 

In this paper we present high resolution (R=20000-45000, cor- 
responding to 14 km/s at 4200A and 6.6 km/s at 9000A) 
spectroscopy of the optical afterglow of Gamma Ray Burst 
GRB050730, observed by UVES@VLT ~ 4 hours after ttie trig- 
ger 

We confirm that the redshift of the host galaxy is z=3. 96764. 
Four intervening systems between z = 3.56 and z = 1.77 have 
been identified along the GRB line of sight, in addition to the 
main system at the redshift of the host galaxy. 

For what concerns the main absorption system, the spectrum 
shows that the ISM of the GRB host galaxy is complex, with at 
least five components contributing to this main absorption sys- 
tem at z- 3.967. Such components are identified in this paper 
with progressive numbers with decreasing velocity values. 

These absorption lines appear in different flavors, from neu- 
tral hydrogen (Ly-a, Ly-/3, Ly-y), to neutral metal-absorption 
lines (O I), and via low ionization lines (C II, Si II, Mg II, Al II, 
Ni II, Fe II, P II, S II) to high-ionization absorption features 
(C IV, SiIV,NV, OVI). 

In addition, we detect strong C II*, Si II*, O I* and Fe II* fine 
structure absorption features in components 2 and 3 (no Fe II* in 
the last component). The ratio of these lines to their ground states 
are used to constrain the gas density and temperature, under the 
assumption that they are produced via collisional excitation, or 
alternatively by UV pumping. If collisional excitation is at work, 
we find that the gas producing the fine structure features has a 
temperature of a few thousand K, and an electron density = 
lO^' + lO* cm""* (component 2) or He = 10+ 10^ cm"-' (component 



Table 4. MetaUicity 
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HI 


22.05 ± 0.29 


- 
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CII 


15.23 ±0.15 








CIV 


15.19 ±0.09 








C 




15.51+0.12 


-6.54 ± 0.41 


-2.93 ± 0.41 


OI 


SAT 








OVI 


SAT 








Si II 


14.96 + 0.15 








Si IV 


15.55 + 0.12 








Si 




15.65 ±0.13 


-6.40 ± 0.42 


-1.94 + 0.42 


Fell 


15.39 ±0.10 


15.39 ±0.10 


-6.66 ± 0.39 


-2.11 ±0.43 



Total column density of the element 



3). On the other hand, if UV pumping is the process responsible 
for the fine structure features, the radiation field exciting the gas 
has an intensity of G/Gq = 10^ -h 10^ 

We speculate that the excitation mechanisms of the fine 
structure lines may be different, with component 3 probably ex- 
cited by collisions, and component 2 by UV pumping. This is 
because the observed O I*/0 1 ratio in component 3 would imply 
a similar Fen*/Fen ratio in the case of UV pumping, which is 
not observed. On the other hand, the radiation field intensities 
calculated from the fine structure ratios for component 2 are aU 
in agreement with each other. 

Component 1, the one with the longest vawelenght, highest 
positive velocity shift, does not present any fine structure or low 
ionization lines; it only shows very high ionization features, such 
as C IV and O VI, suggesting that this component is very close 
to the GRB site. 

For the first three components we can derive information on 
the relative distance to the site of the GRB explosion. Both for 
coUisional excitation and UV pumping, component 2 appears to 
be 10-100 times closer to the GRB site than component 3. 

As to the metallicity, for the system at z=3.967 we obtain 
mean values « 10"^ of the solar metalhcity or less. However, 
this should not be taken as representative of the circumburst 
medium, since the main contribution to the hydrogen column 
density comes from the outer regions of the galaxy while that 
of the other elements presumably comes from the ISM closer to 
the GRB site and since the correction for dust depletion is highly 
uncertain. It is however interesting to note that the mean [C/Fe] 
ratio agrees well with that expected by the single star-formation 
event models of Pipino and Matteucci (2004, 2006). 

We could also measure the [C/Fe] of the single components 2 
and3.The |C/Fe| ratio of component 2 is -0.15+0.13, marginally 
different from that of component 3, 0.54+0.23. Although the re- 
sult is significant at only 2.6cr confidence level, it is interesting to 
speculate on its origin. If the cloud responsible for the absorption 
component 2 is closer than the cloud producing component 3 to 
the GRB site, we expect that dust destruction in this cloud, due 
to the interaction of the GRB radiation field with the medium, 
is more efficient than in the cloud responsible for component 3. 
Since silicates tend to be destroyed more efficiently than graphite 
(Perna, Lazzati & Fiore 2003), we expect more iron and silicon 
freed in the gas phase in cloud 2 than in cloud 3, as observed. 
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